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Abstract 

A surface theoretic view of non-perturbative quantum gravity as "spin- 
foams" was discussed by Baez. A possibility of constructing such a model 
was studied some time ago based on (2+1) dimensional general relativity as a 
reformulation of the Ponzano-Regge model in Riemannian spacetime. In the 
present work, a model based on (3+1) dimensional general relativity in Rie- 
mannian spacetime is presented. The construction is explicit and calculable in 
details. For a physical application, a computation formula for spacetime vol- 
ume density correlations is presented. Remarks for further investigations are 
made. 

1 Introduction 

A surface theoretic view of non-perturbative quantum gravity as "spin-foams" was 
discussed by Baez QTJ] . A possibility of constructing such a model was studied some 
time ago based on (2+1) dimensional general relativity as a reformulation of the 
Ponzano-Regge model in Riemannian spacetime 0. In the present work, a model 
based on (3+1) dimensional general relativity in Riemannian spacetime is presented. 
The main difference is that (3+1) general relativity is not a topological field theory 
but has local degrees of freedom while (2+1) general relativity is identical to the 
non-degenerate sector of BF theory (a topological field theory) in (2+1) dimensions. 
The construction is explicit and calculable in details. This model allows one to in- 
vestigate the difference between the quantum theories for general relativity and BF 
theory from a surface theoretic point of view. For a physical application, a compu- 
tation formula for spacetime volume density correlations is presented. For a family 
of models allowing spin-foam interpretation of quantum gravity, see the references 
HI, H, H, £|, [|, ||. In particular, surface theoretic views are discussed by Baez [p], 
Reisenberger ||, Reisenberger and Rovelli j| and the author In the present 

work, the calculation technologies developed in the connection representation 0] and 
the loop representation M of non-perturbative canonical quantum gravity play im- 
portant roles. 

-'The kind of surfaces playing important roles is called in terms of different names by different 
authors. The name "spin- foam" is due to Baez. 
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Among the models revealing aspects of the spin-foam interpretation, the philoso- 
phy of the present model is similar to that of the Reisenberger model || in the sense 
that both models are based on a path integral of general relativity. The other models 
are based on lattice topological field theories || [|], canonical quantization of general 
relativity Q or others 0. However, the technical difference of the present model 
from that of the Reisenberger model is that the present model uses multiple pairs of 
lattices. In each pair, one lattice is dual to the other. The pairs of lattices are used 
to regularize the exterior product of two forms. The exterior product of two 2-forms, 
which is present in the action functional of general relativity in 4-dimensions, is natu- 
rally defined on the lattice at the intersection of a face and its dual face. In addition, 
the degenerate and non-degenerate sectors of the theory are clearly characterized by 
the number of pairs of lattices. 

In a limit, the partition function of the model is reduced to a product of copies 
of the partition function of lattice BF theory, each of which is independently defined 
on one of the lattices. 2 This fact implies that the contribution from each lattice 
is decoupled in this limit and the use of the multiple lattices for topological field 
theory is redundant. In general, however, the contribution from each lattice cannot 
be decoupled and the partition function of the model is different from that of lattice 
BF theory. This difference is supposed to correspond to the fact that general relativity 
contains local degrees of freedom while BF theory does not. In this way, the model 
provides a way of clarifying local degrees of freedom of quantum gravity. 

The present model is formulated such a way that the inclusion of local degrees 
of freedom of general relativity, which is absent in BF theory, is incorporated to 
coefficients in the partition function and the "trivialization" of the coefficients leads 
to (a product of copies of) lattice BF theory. The computation of the coefficients 
is done by the regularization in terms of the multiple pairs of lattices. Each pair 
consists of two lattices dual to each other. The "trivialization" corresponds to the 
limit mentioned above. 

In the next section, the partition function of the model is presented at the begin- 
ning, and then the rest of the section is devoted to the determination of coefficients 
in the partition function, followed by the presentation of a computation formula for 
spacetime volume density correlations. Remarks for further investigations are made 
in the conclusion section. 

2 The model 

The partition function of the model is 

P „ path 

zpaz\--=u n du( e )x n E G &]xi(n^( e ))> w 

p=leeA p A; J feApA* j eef 

2 The partition function of lattice BF theory in 4-dimensions denned on a single lattice was written 
by Ooguri [§. 
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Cjil[t}:=T,WpA{lhm {j) ({l}). (2) 

Here, (3 is 1 or i (the imaginary unit number) and A is a real or imaginary number. 
£ is an external field coupled to spacetime volume density, / d 4 Xy/g£(x). This field 
is a mathematical tool to calculate spacetime volume density correlations. P is the 
number of pairs of lattices and A p and A* are a pair of lattices dual to each other. 
The index p runs from 1 through P. A family of spin- foams are defined in terms of 
each of the lattices. More details of the lattices and spin-foams are described below. 
Xj(U) is the character of SU(2) element U in the spin-j representation. Y\ path means 
the path ordered product. dU is the Haar mesure for SU(2) element U. e £ A p and 
f E A p mean that they are an edge (1-cell) and a face (2-cell) of A p respectively, 
e £ / means that e is an edge shared by / and other faces; in other words, e is a 
boundary of /. j and I are spins associated with the faces taking values 0, |, 1, 
f • ■ •• W/3,a({^})£) an d are coefficients to be determined below. A particular 

property of these coefficients are as follows. When (5 = 1, A = and £ = 0, then 
£%o[0] = ]lp=i n/GA p A*(2j l + 1) with j associated with each /. Accordingly the 
contributions from the different lattices are decoupled and the partition function 
becomes 

P „ path 

^,o[o] = n n du(e)x n E^'+^di^)) 

p=le£A p A* J feA p A* j ee/ 

2P 

(3) 

_eeA" feA j ee/ 



path 

n jdu(e)x nE(2i+i)xi(n^)) 



Notice that this is not just the partition function of lattice SU(2) BF theory in (3+1) 
dimensions, a topological invariant, but the multiple power of it. We have used the 
fact that the partition function of lattice BF theory does not depend on the choice of 
lattice and let A represent one of the lattices. A line of studies of lattice BF theory 
is detailed in the reference 0. 

The lattices A p and A* used in the partition function are defined as follows. They 
are 4-dimensional piecewise linear cell manifolds consisting of 0-cells (vertices), 1-cells 
(edges), 2-cells (polygons), 3-cells (polyhedrons) and 4-cells. They are dual to each 
other. That is, every k-cell of A p (and A*) intersects with a (4-k)-cell of A* (and A p 
respectively) at a point inside the cells. We call the intersection of a k-cell of A p and 
a (4-k)-cell of A* the "center" of each of the cells. Every cell of A p and A* has one 
and only one center. We define the two lattices such that the center of a k-cell is not 
on the boundary of but inside the cell. The exception is that the center of a 0-cell is 
identified to itself. A spin-foam is defined as a 2-dimensional sub-complex of A p or 
A* with 2-cells labeled by spins and 1-cells labled by intertwiners. A 2-cell labeled by 
the 0-spin is understood as the absence of the cell in the 2-dimensional sub-complex 
under consideration. 

In the model, we use the multiple of such pairs of lattices and denote the set of 
the lattices by {AA*} meaning {A p , A*{p — 1, • • • , P)}. In addition, we impose two 
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conditions to the lattices. One condition is that the position of the center of every face 
of A p coincides with the position of the center of a face of A q for all p, q = 1, • • • , P. 
This condition is needed in order to characterize the degenerate and non-degenerate 
sectors of the theory. The other condition is that A p and A* are isomorphic to A q and 
A* respectively for all p, q — 1, • • • , P and hence the faces of them correspond one-to- 
one. In other words, A 1; • • ■ , Ap are copies of the same (abstract) lattice embedded 
in different locations of the spacetime manifold with a restriction due to the first 
condition and so are their respective dual lattices A£, • • - , Ap. The positions of the 
one-to-one corresponding faces of A p and A q under the isomorphism do not necessarily 
coincide with each other. The second condition regularizes the kind of lattices utilized 
for the model and guaratees that copies of lattice BF theory independently defined 
on single lattices produce an identical partition function. (A way of constructing the 
pairs of lattices is discussed in Appendix [A].) 

In order to determine the coefficients Wp^({l},^) and ^'({l}) we compute a 
path integral of general relativity. The action for general relativity we consider is 
the Plebanski action [|10| . It consists of two terms, both of which is spacetime dif- 



feomorphism invariant separately in addition to their internal gauge invariance, and 
is 

S[A, B, 4>] := J d'x^r hcd B\ h Fi d - ±j d*x^~e abcd B ab Bi d . (4) 

Here a, b, c ■ ■ ■ are spacetime indices {0, 1, 2, 3} and i,j--- are internal SU(2) adjoint 
indices {1,2,3}. e abcd is the alternating tensor with density weight 1. F^ b := d a A\ — 
dbA a + e l j k A J a A b are the curvature of an SU(2) gauge connection A a , B ab is an SU(2) 
algebra valued 2-form field and 4>ij are scalar fields (traceless symmetric with respect 
to the SU(2) indices). The dimensions of A l a , B l ab and </>y are of length inverse, length 
inverse squared and length inverse squared respectively. A is a parameter with the 
dimension of length squared and its presence can be understood if one rescales A l a 
and B ab appropriately. The Euclidean sector is obtained by real valued fields and 
imaginary time coordinate. The five scalar fields 0^ are Lagrange multipliers forcing 
the non-trace parts of e abcd B l ab B 3 cd (about the internal indices i and j) vanish. This 
condition implies the exsistence of non-degenerate tetrad fields e° a and e l a such that 
B ab = e^e b + \e % jk e J a e b if 5ij€ abcd B l ab B J cd ^ 0. In other words, in its non-degenerate 
sector the action has the value of the self dual Hilbert- Palatini action evaluated on the 
tetrad. The first term of (|j) is known to be the action of the BF theory, a topological 
field theory. The second term amounts to introduce local degrees of freedom contained 
in general relativity. 

We define a path integral in terms of an action projected to the pairs of lattices 
{AA*}. The variables projected to the lattices are defined as follows. For an edge 
e, define U(e) := joexp[/ e dx a A a T i ] 1 a path ordered parallel transport along e. Here Ti 
is the Pauli matrix divided by 2i. For a face /, define rf{f) := | dx a dx b B l ab . The 
scalar fields (fiy are defined at the centers of faces. The integrals over edge and face 
are oriented such that they are compatible with the orientation of A„ and A* 
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The action projected to the pairs of lattices {AA*} is defined by 

P I i path 

S {AA * } [U, V , </>}:= £ E [vVhnriU(r)I[U(e)U-\n} 

P=i/eA p A; V e e/* 

~Mfn[vV)rf(n+rf(f)vV*)^ , (5) 

where f* is the face sharing its center with / on a pair of lattices A p and A*, and //* 
means the center of / and f*. U(f*) is a parallel transport defined not on an edge 
of the lattice but on a curve from the center of /* to one of the vertices belonging 
to /*. The choice of the vertex is arbitrary and does not affect the definition of the 
partition function of the model. This parellel transport is necessary in order to make 
the first term gauge invariant since the two faces / and /* intersect with each other 
only at their center and the untraced parallel transport must base on the center. 

The exponential of the action can be considered as (an extension of) the graph- 
cylindrical function discussed in 0. The integral measure we use for the variable A l a 
is the Ashtekar-Lewandowski measure f7|. The measure for the scalar fields <pij is the 
one discussed in ||. The measure for the variable B l ab is analogously defined. The 
path integral has the form of 



dti(A)d(i(B)dn(<j>)V { AA»},ii,( A > B > 0) := 

/ n d u(e) n d v(f) n w/> 

ee{AA*} /e{AA*} /e{A} 

V(*7(ei, A), ■ ■ ■ U(e h A); V (f h B), ■ ■ ■ v (f m , B); (j>(vi), ■ ■ ■ cj>(v n )), (6) 

where ^{aa*},^ is a cylindrical function defined on the pairs of the lattices {AA*} in 
terms of a complex valued integrable function ijj on [SU(2)] 1 x [sw(2)] m x R n . Note 
that e G {AA*} does not mean that e is a lattice but means that e is an edge of one 
of the lattices. In the same way, / G {AA*} means that / is a face of one of the 
lattices. We loosely use this kind of notation. This integral form is diffeomorphism 
invariant. 

The path integral which will be identified to (HD is 

z'eM ■= fdu r d ^[^->^+^E/eA,wr)*^(r)] 

J J —oo 

/roo 
dU(e)x n / dr,(f)x 
eetAA-J- fe{AA*} J -° D 

/OO 
#(//*) 5{<t> 1X + <P22 + 033) X 

Tew - 00 

n n /^(f)e* (/1 ^ W)I ^ [ ' (erl(/,)1 x 
P =i f£A p A; 

P 

rr yi e -¥ x ^(ff*)W(f)riHfn+riHfW(fn] e mfrHU)vHn_ ^ 

p=l /eA p 
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Here, the external field £ is defined at the centers of faces as 4>ij is defined. We note 
that if P = 1 then one can easily solve the constraints on rf imposed by (fiij and finds 
only degenerate solutions corresponding to 5iji abcd B l ab B J cd = 0. In order to include the 
non-degenerate sector, one needs P > 1. (See Appendix [B] for more discussion on the 
non-degenerate sector.) 

Let us expand the exponential containing the trace of an SU(2) element to its 
characters as follows. 

J dVe "J w W~ x 

2 3 (-l) h+h+h (2h + 1)(2Z 2 + l)(2/ 3 + 1) x 

Xh invuv- 1 ]^ [tzVUV- 1 ]^ {nvuv- 1 } 

J2h+i(/V) J2i 2 +i{f3rf) J 2 z 3 +i(/V) „ ^ 
= ,S, _ 8? W «-"ta( f/ >' < 8 > 

Here, J m (x) is the Bessel function of the first kind and we have used the following 
formulae. 

e^™ = ^2^±±iVj 2j+1 (x) Xj (U), (9) 

3 X 

2lT J-1T 

The first formula can be easily proved and the second is a definition of the Bessel 
function. ^'(h,h,h) is defined as follows. 

n {j) (h,l 2 , k) := 2 3 (-l) il+ ^(2/ 1 + 1)(2Z 2 + l)(2/ 3 + 1) x 

E ^(l)^i 1 (n)^ 2 (r 2 )^i 3 (r 3 )(-l)" 3 -" 13 E (2iu + l)x 

mi,rii,m,n Ji2,mi2,ni2 

/ Zl Zl2 \ Ml Z 2 Zia \ / Zl2 ^3 j\ / *12 ^3 J \ 

Vm n 2 ni2/ VTOi m 2 m i2 / \ — «i2 n 3 n)\-m 12 m 3 m ) ' 



Here, D^J n (U) is the spin-j representation matrix of SU(2) element U and m and 

n run from — j through j with the increment 1. ( 31 ^ 2 ) is the so-called 

V mi m 2 m J 

3j-coefficient and defined by 

( 31 32 3 ):= L}fL^Jl^ imi . j 2m2 \j m ) (12) 

with the Clebsch-Gordan coefficient (jimi; j' 2 ot 2 |j'to) and we have used properties of 
D^JJJ) such as 

(-l) n - m D { £(U) = D U) n % m (U), (13) 
J dUD£l(U)D2:,(U) = ^^<^ W <W, (15) 



where the asterisks mean the complex conjugate and the sum of j is taken over 
\ji ~ J2I through ji + j'2 and the sums of n and m over — j through j. Then the 
partition function can be written as follows. 



//*oo 
dU(e)x ]J / Mf)x 

efc^ZA-j- /e{AA«}^°° 

/oo 
#(//*) 5(0n + 22 + 33 ) x 
-00 



X 



/e{A} 

yr tt V 4(/'Hi(^(/)) %') + i(M/)) JwiW(/)) x 
M/el{%/, W) 

W(D) J 2W)+ i(/W *)) 

/W*) /W*) /W*) 

e - fWij (//* ) fo* ( W (/* ) HV (/)»,* (/* )] e m(ff*)Vi (fW (/* ) x 

path 

3(f) eef 

path 

E n (, ' (/,)) (ii(r),i2(r),/s(r)) W . ) (ii ( 16 ) 

Here, {/}//* stands for h(f),l2(f),h(f),h(f*),k(f*),h(f*)- Let us write the other 
exponentials explicitly and expand them as a power series as follows. 

e - 1 /JAtfy (//* ) to* W (/* (/)»?* (/* )] g W/* )th (fW (/* ) 

= e -W 11 (//*)r ) 1 (/) f | 1 (/*) e -W 2 2(// , )/) 2 (/)r/ 2 (/*) x 

e -W 3 3(// , )^(/)0 3 (/') e -Wi2(//*)[» 1 (/)'/ 2 (/ , )+') 2 (/)') 1 (/*)] x 

e -i/3A0 23 (//* ) [»? 2 (/)r, 3 (/* )+r? 3 (/)r? 2 (/* )] g-i/3A</> 3 i (//* ) [v A (f)v 1 (f*)W (f)v 3 (/* )] x 

e Wr ) fo 1 (/h 1 (/* )+»j 2 (f)v 2 (/* )W (fW (/* )] 



N u \ Mi! 



E E 1 " TV n WifWin 



1/ r*\lJVii+Mi 



X 



ATn=0Mi=0 

.V, ,! ~" ' Mo 



E E 1 ;J r "7V J fr 2 (/V(/*)r 22+M2 

7V 22 =0M 2 =0 iV 22 

Ar 3 3=0M 3 =0 iV 33- lvl 3- 

E E [ w £ ,1/ [^(/Wn^VWcrr 21 x 

jV 12 =0JV 21 =0 iVi 2 !iV 2 i! 

7V 23 =0Ar 32 =0 iV 23 !7V 32 ! 

E E [-^({^p 1+JVl3 [ r/ 3( / ) r/ l( r) ]^ [7/ l (/)?7 3 (r)ri 3. 

AT 31 =0Ari3=0 JV31.7V13. 



X 
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Here, the non-negative integers A^- and with i,j = 1, 2, 3 are associated with the 
pair of faces / and /*. Then the partition function becomes 

Z'fiM\ = II [dU(e)x{[ I] E %(W//«.Ox 

eG{AA*} J p=l/eA p {(} /r 

path 

T,^ if Hh(f),i2(f),h(f))xm(I[u(e))x 

j(f) eef 

path 

with 



Wp,x({l}fM) ■= / Mf)dv(r) / #(//*) 5(0n + 22 + 033) x 

J —oo J — oo 

HWf [-w 22 (//*)]^ [-w 33 (//1* 

{Af,M} /; * 



^ A;,! iV 22 ! .V,,! 



x 



HMjr; HWf!^ [-j^(/rF a+Jfu 

A^ 12 !iV 21 ! A^Ay Ar 31 !7V 13 ! 

[^(//*)] Ml+M2+M3 
.\/;!.\/,!.U ;! ! X 

|-^l^^jMi+JVii+JVi2+JVi3 ^2^^M 2 +Af 2 2+A r 2i+A r 23 ^3^jM 3 +JV33+JV 3 2+.N3i x 

^l^*^Mi+jVil+AT 21 +jV 31 ^2^*^jM 2 +JV22+JVl2+JV32|-^3^*^M3+JV33+JV23+JVi3 x 

4(f )+ iW(/)) <W)+i(/V(/)) 

/W) /W) /W) 
j2t 1(/ ) + i W(D) ^(fl+iW)) ^(/hi 

W) /W*) /W*) 



X 



(19) 



Here, {A^}//* stands for the non-negative integers iV^ and Mj associated with the pair 
of faces / and /*. The integrations over rf can be performed. Here, for concreteness, 
we compute the coefficient for the case that (3 is positive real number. A useful 
formula is J(m,n) := $ cx ' 00 x n ~ l J m {x)dx = (3 n Jf^ x n ~ x J m {(5x)dx and 



J(m, n > 1) = i 



n-l 



d 1 \ p~ im6 \ + ( — - ^ g-imSl 

d9 cos 6 16=0 \d9 cos 9 



,(20) 



J(m,0) = -[l-(-in. (21) 



Here 4z — a means 4z — 5 acts n — 1 times. For given m and n, Jim, n) exists 

\ (19 cos 9 J d9 cos 9 & ' V ' / 

and can be explicitly computed. In terms of J(m,n), the coefficient Wp t \({l} ff*, £) 
can be written as follows. 



/oo 
#(//*) ^(011 + 022 + 033) 
-oo 



X 
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X 



[-ip-^Mff*)]" 11 HElMdll¥l [-iP- 1 *Mfn] Nm 

{ »M} tr ^ M 

[-ip-'xMr)]" 12 ™ 1 H/r^M//*)]" 23 *" 32 Hgjg//t^ 

iV 12 !iV 21 ! 
[^ _1 e(//*)] Ml+M2+M3 



N 31 \N 13 \ 



x 



Mi!M 2 !M 3 ! 

J(2Zx(r) + 1, M 1 + jV u + iV 12 + N 13 )J(2l 2 (f) + 1, M 2 + iV 22 + jV 21 + iV 23 ) x 
J(2Z 3 (/*) + 1, M 3 + iV 3 3 + iV 32 + N 31 )J(2h(f) + 1, Mi + iV n + iV 21 + iV 31 ) x 
J(2/ 2 (/) + 1, M 2 + iV 22 + AT 12 + N 32 )J(2l 3 (f) + 1, M 3 + iV 33 + 7V 23 + iV 13 ). (22) 

Now, define the coefficients and Wp,\({l}-,£) such that Z'p t \[£\ is identi- 

fied to ^/3,a[C] as follows. 

n w ({i}):=n II 2-*n U{f)) {h{f)Mf)Mf)), (23) 

P =l /eApA* 

^,a(w,o :=nn ^, A ({/} /r ,o. (24) 

p=l /6A P 

Here, 2~ 9 and /3 6 have been inserted for convenient normalization. (|23|) with (|TT|), 
and (^4]) with (fffi) are the final forms for the coefficients (§) in the partition function 
of the model. 

For a completeness, let us compute the coefficient Ci^[£] for the case that (3=1, 
A = and £ = 0. In that case, 

VMW,0) = f[ II J(2/i(/) + l,0) J(2Z 2 (/) + l,0) J(2Z 3 (/) + l,0) 
p=i /ga p a; 



p 

n n 

p=i/sA p a; L 



■1) 



2/i(/)+l 



2Zi(/) + l 



1 _ (_l)2«2(/)+l 

Mf) + 1 



1 _ (_l)2J 3 (/)+l 

2M7) + 1 



(25) 



and 



cS[o] = E^i,o(W,o)^ ) ({/}) 

{0 

= n n Et 1 - (-if - (-i) 2 ' 2+i ][i - (-i) 2 '3+i] 2 - 6 (-i)'i+^3 

p=i/eA p A; h,h,i 3 

E ^(l)^l(n)^i 2 (r 2 )^i 3 (r 3 )(-l)^- m3 E (2Zi2 + l)x 

'l2,"ll2,ni2 



X 



mi,ni,m,n 





z 2 






n 2 


n 12 J V 



Zo z 



12 



kz h j 
-n 12 n 3 n 



Z12 Z 3 j 
-m 12 m 3 m 



= 11 II E Xh(n)xh(T 2 )xi 3 (T 3 )xi 12 (i) x 

p=i/eA p A; h,h,hM 
J dVdW Xh {nV) Xh (r 2 V) Xh2 {VW-^xis {rsW)xj{W) • 



(26) 



Here, we have used the formulae that Xiiji) — \{~ 1)' if I is an integer and Xi( T i) = 
if / is a half- integer and that =21 + 1 and that 





32 




r 


n 2 





Jl ]2 33 
mi 771,2 ^3 



dVD^{V)D%l 2 {V)D%l 3 {V). (27) 



Let us work out the coefficient C^o[0]. From the last expression, V and W are forced 
to be V = W = 1 after the sums over 1%, l 2 , l 3 and l\ 2 and the integrations. Then we 
get 

csn°] = n n ^(i) = n n ( 2§ ) 

p=i /gApA* p=i /gA p a; 

if j is an integer; otherwise, vanishes. Here, j is associated with each /. This "trivial" 
coefficient makes the partition function of the model the multiple power of that of 
the lattice BF theory as claimed already. 

Finally, for a physical application, we present a computation formula for space- 
time volume density correlations. From the partition function of the model, define 

zfx" Xn) m which is Zp t \ [£] with w MU) re P laced b y sfo-'-^wWtfU)- 

Then the spacetime volume density correlation function is 

(V^"-v£M) ; = Z ^ rni [ ° ] - ^ 

This quantity is the average of a product of spacetime volume densities defined at the 
centers of faces with respect to the partition function. The spacetime volume density 
at x is computed by 



g(x) := £ ^(/M/*), (30) 

ff*=x 

where the sum is taken over faces / and /* whose center //* is located at x. 



3 Conclusion 

We have presented a model for quantum gravity, which could reveal an aspect of the 
spin-foam interpretation of spacetime. Its partition function is defined by ([[]) with 
coefficients @. We have determined the coefficients based on a path integral of (3+1) 
dimensional general relativity in Riemannian spacetime. To do so, we projected the 
Plebanski action (|J) to multiple pairs of lattices, each pair of which consists of two 
lattices dual to each other, (Bp and then computed its path integral (^). The final 
forms of the coefficients are given by ( p3|) with (fOP , and ([24]) with (0). In terms 
of the partition function, we have presented a computation formula for spacetime 
volume density correlations. 

We conclude with remarks for further investigations: 
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(i) The variables U and </>y are left unintegrated. The integration for U is known 
to produce functions of spins and intertwiners labeling respectively the faces and 
edges of a spin-foam. This integration is already present in the case of the BF theory 
and the resulting functions consist of the so-called 3nj-symboles. A spin-foam can be 
decomposed to a family of surfaces and these functions can be understood as telling 
the number of "distinguishable" ways with respect to a given lattice (dual to the 
lattice the spin- form bases on) the family of surfaces are allocated as studied in ||. 
The integrations for (f)^ are particular in quantization of general relativity and are 
supposed to induce local degrees of freedom. In the present model, physics of the 
local degrees of freedom may be translated from the information how single spin- 
foams "interact" with each other. The coefficient Wp t \({l},£) has roughly the form 

of J dxi ■ ■ ■ dxk J2 ni -n k fni-n k with x' s corresponding to 0^ and naive term- 

by-term integration leads to divergences. We have to understand the properties of 
fm-n k before the integrations of 4>ij- 

(ii) To perform the path integral of general relativity, we have used diffeomorphism 
invariant measures. However, we have to check the resulting partition function is 
really diffeomorphism invariant. In a limiting case, we have shown that the partition 
function is reduced to the multiple power of the partition function of the lattice BF 
theory, which is a topological invariant and certainly diffeomorphism invariant. 

(iii) A technical peculiarity in the construction of the present model is the use 
of multiple pairs of lattices. Each pair consists of two lattices dual to each other. 
Because of this, physics comes out not from single spin-foams but from "interactions" 
between them. We would like to understand whether the use of the multiple lattices 
is necessary or redundant. 

(iv) In the standard lattice gauge theory, in which the number of degrees of freedom 
is finite and each variable is compactified, the use of lattice reguralizes the divergences 
due to the integrals along the gauge orbits as long as computations of gauge invariant 
quantities are concerned. However, in the present model, the reguralization is not yet 
under control because not all the variables are compactified. In particular, to perform 
computations of gauge dependent quantities such as propagaters, the knowledge of 
how to fix a gauge is crucial. 

(v) As physical applications, one wants to compute physical quantities. For ex- 
ample, a detailed study of the spacetime volume density correlation function we have 
presented is to be pursued. 

(vi) In order to understand the physical meaning of spin-foams, a possible way of 
studying it is the following. Couple external fields to the variables A l a , B l ab and 4>%j and 
then compute the partition function with the external field. In terms of the partition 
function, one may define Schwinger equations corresponding to the classical Einstein 
equations with those variables. Detailed study may allow one to understand the 
physical meaning of the spin-foams corresponding to classical spacetime geometries. 

(vii) Ultimately one wants to formulate a quantum gravity model in Lorentizn 
spacetime. As an extension of the present model, Lorentzian spacetime gravity is 
defined in terms of the path integral of complex variables with certain reality condi- 
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tions. A possibility of including the reality conditions into the path integral ought to 
be studied. 

A Pairs of lattices 

We discuss in this appendix a way of constructing the multiple pairs of lattices de- 
scribed in Sec.^. Because 4- dimensional lattices are difficult to visualize, we first 
utilize 2-dimensional lattices , the role of whose edges and dual edges mimics the role 
of faces and dual faces of 4-dimensional lattices, to construct pairs of lattices. Then 
we repeat the construction in 4-dimensional case. Note that these constructions are 
not oriented to show a mathematical proof of the existence of the pairs of lattices. 

On a 2-dimensional manifold with the topology of S 2 for simplicity, draw a lattice 
and its dual lattice. They consist of 0-cells (vertices), 1-cells (edges), and 2-cells 
(polygons or faces). The edges are not necessarily straight lines but may be curved. 
Every edge of the lattice (and the dual lattice) intersects an edge of the dual lattice 
(and the lattice respectively) at their "center." Let us plot a dot at every center 
of edges (and dual edges). Draw curves connecting the every two centers of edges 
sharing a vertx of the lattices. This drawing creates another lattice, whose vertices are 
located at the centers of the edges of the original lattices. We call the created lattice 
"child" of the original lattices. Erase the original lattices. We have now only the 
child on the 2-dimensional manifold with dots indicating the positions of the vertices 
of the child. 

Move the child as follows while the dots are fixed. Every vertex is moved from one 
dot to another or is not moved so that every dot is occupied by only one vertex. The 
edges sharing a vertex are moved together with the vertex if the vertex is moved. Let 
us call the resulting lattice "cousin" of the child. They are isomorphic to each other. 

Draw a pair of lattices (a lattice and its dual lattice) from the cousin such that 
their relationship is analogous to the relationship between the original pair of lattices 
and their child. Erase the cousin. Now, we have only the second pair of lattices with 
dots on the 2-dimensional manifold. The centers of edges of this pair of lattices are 
located at the dots and hence the position of the center of every edge of this pair 
of lattices coincides with the position of the center of an edge of the original pair of 
lattices. One of the original pair is isomorphic to one of the second pair while the 
other lattice of the original pair is isomorphic to the other lattice of the second pair. 
By repeating the same process one can generate a third pair, a fourth pair, and so 
forth. 

Next, let us repeat the analogous construction for 4-dimensional case. In a 4- 
dimensional spacetime manifold with the topology of S" 4 for simplicity, draw a lattice 
and its dual lattice. They consist of 0-cells (vertices), 1-cells (edges), 2-cells (polygons 
or faces), 3-cells (polyhedrons) and 4-cells. The edges and polygons are not necessarily 
straight lines or flat surfaces respectively but may be curved ones. Every face of the 
lattice (and the dual lattice) intersects a face of the dual lattice (and the lattice 
respectively) at their "center." Let us plot a dot at every center of faces (and dual 
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faces). Draw curves connecting the every two centers of faces sharing a vertex of 
the lattices. The curves connecting the centers of faces sharing a vertex form a 
polygon. This drawing creates another lattice. Its vertices are located at the centers 
of the faces of the original lattices. We call the created lattice "child" of the original 
lattices. Erase the original lattices. We have now only the child in the spacetime 
manifold with dots indicating the positions of the vertices of the child. 

Move the child as follows while the dots are fixed. Every vertex is moved from one 
dot to another or is not moved so that every dot is occupied by only one vertex. The 
edges sharing a vertex are moved together with the vertex if the vertex is moved. Let 
us call the resulting lattice "cousin" of the child. They are isomorphic to each other. 

Draw a pair of lattices (a lattice and its dual lattice) from the cousin such that 
their relationship is analogous to the relationship between the original pair of lattices 
and their child. Erase the cousin. Now, we have only the second pair of lattices with 
dots in the spacetime manifold. The centers of faces of this pair of lattices are located 
at the dots and hence the position of the center of every face of this pair of lattices 
coincides with the position of the center of a face of the original pair of lattices. One 
of the original pair is isomorphic to one of the second pair while the other lattice of 
the original pair is isomorphic to the other lattice of the second pair. By repeating 
the same process one can generate a third pair, a fourth pair, and so forth. 

B Non-degenerate sector 

Pairs of lattices in 4-dimensions are used in the model. Each pair consists of two 
lattices dual to each other. That is, every k-cell of one lattice intersects with a (4-k)- 
cell of the other lattice at a point inside the cells. The intersection is called "center." 
In particular, every face (2-cell) of one lattice shares a center with a face of the other 
lattice. These two faces are said to be dual to each other. Let A p and A* denote a 
pair of lattices with p = 1, • • • , P. The P is the number of pairs. If they form the 
multiple pairs described in Sec.0, then there exist a face f p in A p and its dual face 
f* in A* with their center denoted by f p f* such that fif* = ■ ■ ■ = f P f P . (Here the 
equality means that their positions coincide.) 

Let //* represent the position of the coincident centers. The Lagrange multipliers 
4>ij are defined at //*. The constraints imposed by the multipliers at //* are 

X>um/;) + WpM/;)] = c(fns ij . m 
p=i 

Here C(ff*) is an arbitrary non-negative real number determined by rf(f p ) and rj l (f*) 
and proportional to spacetime volume density at //*. C(ff*) > and C(ff*) = 
mean that the solutions for rf are non-degenerate and degenerate respectively at //*. 

It is easy to check that if P = 1, then solutions for rf at //* do not exist unless 
C(ff*) = 0. This fact means that the model needs more than a pair of lattices in 
order to contain the non-degenerate sector. Once the number of pairs, P, is fixed, 
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then the non-degenerate and degenerate sectors can be clearly characterized. Since 
the constraints have the form of the inner product for three vectors (i,j = 1,2,3) 
with their norms proportional to spacetime volume density, the non-degenerate sector 
of the model is analogous to the space spanned by non-degenerate orthogonal vectors 
in a 3-dimensional linear space. 
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